Kondo effect in multielectron quantum dots at high magnetic fields. 
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We present a general description of low temperature transport through a quantum dot with any 
number of electrons at filling factor 1 < v < 2. We provide a general description of a novel Kondo 
effect which is turned on by application of an appropriate magnetic field. The spin-flip scattering of 
carriers by the quantum dot only involves two states of the scatterer which may have a large spin. 
This process is described by spin-flip Hubbard operators, which change the angular momentum, 
leading to a Kondo Hamiltonian. We obtain antiferromagnetic exchange couplings depending on 
tunneling amplitudes and correlation effects. Since Kondo temperature has an exponential depen- 
dence on exchange couplings, quantitative variations of the parameters in different regimes have 
important experimental consequences. In particular, we discuss the chess board aspect of the exper- 
imental conductance when represented in a grey scale as a function of both the magnetic field and 
the gate potential affecting the quantum dot. 

PACS numbers: 72.15.Qm 



I. INTRODUCTION 

Kondo effect has been recently observed in quantum 
dots (QD) coupled to leads by analyzing the tempera- 
ture dependence of the conductanceErtl Kondo physics 
has been mainly detected for an odd number of electrons 
N in the QD while for even N, the temperature depen- 
dence of the conductance is usually that corresponding to 
a normal system. This is interpreted in terms of the for- 
mation of singlets by pairs of electrons so that only when 
N is odd a single electron remains unpaired. This last 
electron, in a state doubly degenerate by spin, is reSjDQii- 
sible for the Kondo behavior propoped-iong time agcoEi 
However, some recent experimentsliinL^I suggests a more 
complicated situation because Kondo effect is also ob- 
served for even N in some situations. A particularly in- 
teresting experimental feature is the appearance, for a 
given N, of alternating high and low conductance valleys 
as a function of an external magnetic field HIJO. When 
N is varied in ±1, the high and low conductance valleys 
are interchanged. Therefore, the representation of the 
conductance (in a grey scale) as a function of both B 
and a gate potential W-hidb allows to vary N, takes the 
aspect of a c/ies ( s_.6parcOllj. 

Three modelsLJ E£l have been already proposed for un- 
derstanding how Kondo effect might be possible for states 
of two electrons. In these approaches, the authors con- 
sider the situation of double degeneracy that a magnetic 
field can create between the singlet (S z — 0) and the 
triplet (S z — 1) states of the electronic pair. However, as 
shown below, the actual situation is, in general, rather 
more complicated due to correlation and tunneling am- 
plitude effects. 

In this paper we present a general description of the 
Kondo effect in a QD in the presence of a high magnetic 
field such that the filling factor is 1 < v < 2. This is-a 
regime in which experiments clearly show Kondo effcctLil. 
The description is valid for any number of electrons (even 



or odd without restriction to 1 or 2) and any value of 
the QD spin S z (not restricted to 0, 1 for even N or 
1/2 for odd N). Instead of describing spin-flip scattering 
in terms of spin-ladder operators we find a set of 

spin-flip Hubbard operators describing a collective spin 
effect of all the N electrons contained in the QD. De- 
spite both N and S z can be very large, the scattering 
of carriers only produces transitions between two many- 
body states of the QD with S z differing in ±1. These 
many-body states are well characterized theoretically in 
the regime 1 < v < 2, and we find that the spin- flip pro- 
cess can be described by a Kondo Hamiltonian with anti- 
ferromagnetic couplings which depend on both tunneling 
amplitudes and correlation effects. As a consequence, the 
system presents a Kondo behavior with a Kondo temper- 
ature Tk which we analyze in various limiting cases. 

In section II we discuss the ground state (GS) of an 
isolated QD in the regime 1 < v < 2. in Section III an 
effective Kondo Hamiltonian is obtained, by means of a 
scattering description of tunneling through the QD. Sec- 
tion IV contains a discussion on the exchange couplings 
and the Kondo temperature in the regime 1 < v < 2. 
Section V is devoted to the explanation of the chess board 
aspect of the experimental conductance. A summary is 
given in section VI. 



II. QD SPECTRUM 

We consider N electrons in the presence of a magnetic 
field and confined in a QD coupled to leads. The Hamil- 
tonian is 



H = H, 



QD 



H T . + H 



TUN- 



(1) 



Within a lowest Landau level approach, an isolated 
parabolic QD with magnetic field and interaction be- 
tween the electrons is described by (hereafter we take 
h = 1) 
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H QD= — + " Y 1 — + 9VbBS z + 

2 ^ t Vm 1 m 2 m 3 m A d\ ni ai d\ n2 a2 d m3 ^ z d m/L ^ i . (2) 

The first two terms describe a single particle contribu- 
tion depending on both the QD confinement loq and cy- 
clotron uj c frequencies through £1 = y/u>% + AuJq . M is 
the third component of the total angular momentum. 
The third term is the Zeeman energy, depending on the 
Lande g-factor. In the last term, describing electron- 
electron repulsion, dj^ CT creates an electron with angular 
momentum m and spin a in the QD. The Coulomb in- 
teraction matrix elements V mi m 2 m 3 m 4 have a typical en- 
ergy scale e 2 /elB, where e is the dielectric constant and 
Ib = l/VmCl the magnetic length. These interactions 
decrease with the increasing width w of the quantum 
well in which the QD has been built up. 

Hl = J2k,(T £ k,<y c k.cr c k,o- describes, in a single-particle 
approach, the leads having electrons with quantum num- 
bers a occupying states up to the Fermi energy ej?. 

A crucial role is going to be played by the tunneling 
part of the Hamiltonian 

Htun = Vm ( d ln^ c k,a + cLfW) (3) 

k.m,a 

in which we neglect any dependence on k of the tunneling 
amplitudes V m (taken as real positive) but we retain the 
dependence on m because is going to produce physical 
consequences. 

Hqd can be numerically diagonalized for a signifi- 
cantly troad range of TV, as discussed in many theoretical 
paperstZI. An example (similar results are obtained for an 
odd number of electrons) is shown in Fig. [I] and ||, fee 
N = 8, which is the regime in which some experimentstil 
have been performed, although some otherscj involve a 
larger number of electrons. In order to have access to all 
the possibilities for the GS, it is convenient to discuss a 
case with strong interaction effects, so that we present re- 
sults for a QD with w — 0. The case of an arbitrary finite 
width w can be qualitatively discussed from these results: 
correlated GS tend to disappear when w increases. Fig. 
[I] shows the phase diagram of the possible GS between 
v = 1 and v = 2. The wider regions of the upper part of 
the phase diagram correspond to compact states which 
are favored by large Zeeman coupling 

A'-l N-K-l 

\c§- K )=u d i,i n 4, T i°> w 

m— m— 

where |0) is the vacuum state. In going from left to right, 
one finds successively the states \Cf) [y = 2), |C|), |C|), 
\C]) and |Cg) (y — 1). Strong interaction effects are 
manifested in the appearance, lower in the phase dia- 
gram, of the narrowiec^gions corresponding to different 
skyrmion-like statest2lli*2 



\SK% tKt± ) =Mp,± (A+ ±1 ) \C§_ k ) (5) 

of topological charge 1. A/p,± is a normalization constant 
and 

AL=D m+1 ) Tl/2 4+i,i<W' (6) 

m 

Adjacent regions correspond to values of P which differ 
in one unit and, as a consequence, their spins also differ 
in 1. Dashed lines in Fig. |l| depict, for two different val- 
ues of wo, the evolution of the GS for a GaAs QD within 
the range of magnetic fields (perpendicular to the QD) 
given at the edges of the lines. 



GaAs QD; N=8 




FIG. 1. Phase diagram of the possible GS of a GaAs QD 
with 8 electrons at 2 > v > 1. The axes contain kinetic 
(fi — uj c )/(2e 2 /eIb) and Zeeman g* = gfisB /(e 2 /eIb) con- 
tributions to the energy. The (upper) wider regions corre- 
spond to compact states from |C|) (y = 2) at the left to \C%) 
(y = 1) at the right. The (lower) narrower regions correspond 
to different skyrmion-like states. Dashed lines depict, for two 
different values of the QD confinement frequency u>o, the evo- 
lution of the GS within the range of magnetic fields given at 
the edges of the lines. 

Fig. H shows the evolution of the GS properties 
from v = 2 to. v = 1 for the system of Fig. [j] with 
c^o = 5.4meVtJ. The GS energy Eqs has a kink any 
time a crossing of states occurs. Since these kinks are 
not obvious in Fig.|^(b), we represent, in part (a), the 
energy splitting AE = E exc — Eqs between the lowest 
excited state and the GS. The spin S z and the third com- 
ponent M of the total angular momentum of the GS are 
also given in (c) and (d) respectively. 
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GaAs QD; co n =5.4meV: N=8 
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FIG. 2. Evolution with the magnetic field B of different 
magnitudes for the GS of a GaAs QD with TV = 8 and 
luq — 5.4meV. Energies are measured in units of e 2 /eIb- The 
left side of the figure corresponds to v = 2 and the right side 
to v — 1. 

Independently of the particular details, there is a gen- 
eral trend for any N: when the magnetic field increases 
from v = 2 to v = 1, the GS changes many times 
(at least N/2 changes corresponding to electrons N/2 
flipping their spins one by one). In other words, for 
2 > v > 1 there are several crossings of the lowest energy 
states. This means that, for many values of the field, the 
GS is doubly degenerated (or non degenerated but having 
an excited state with extremely low excitation energy). 



Moreover, in almost all the cases, these two states have 
spins differing in 1 while they have different M. 



III. KONDO HAMILTONIAN 

As a consequence of the above numerical result, the 
main physics of the problem is captured by a two level 
system approachEil. One considers the QD in the regime 
2 > v > 1 as having two degenerate (or almost degener- 
ate) GS's \GS?) and \GS¥) with M t ^ My, and spins 



differing in 1, i.e. (GS$\S Z \GS$) = (GS^\S Z \GS?) + 1 
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in which 



Htun mixes these GS's with states \N ± 1, 
N ± 1 electrons are in the QD. In our two level system 
description, the tunneling Hamiltonian is projected on 
the subspaces subtended by the two GS's for N elec- 
trons and the connecting |JV± 1) electron states. In this 
process, we consider the connecting \N ± 1) states as 
non-degenerateE3. Using the notation S = {tf,-IJ-}! we 
introduce the tunneling spectral amplitudes 



A_, E = XV OT <JV-l|d m , ff |GS : 

in, a 

A +tI1 = J2Vm(N + l\dljGS : 

m,<7 

and the spin-flip Hubbard operators 
Xv t v = \GS»)(GSg\. 



(7) 



(8) 



From projected tunneling Hamiltonian Htun, the inter- 
action between the QD with N electrons ancL-the leads 
is studied by a standard scattering description^. An ef- 
fective coupling is obtained by summing up to all the 
possible intermediate states |7), in a second order per- 
turbation approach: 



H, 



eff 



E 



E 

fc.fc'.E 



Htun\1){I\Htun 

E GS - E I 
A 2 s 5fe,fe'Xs,E 



A 2 r 



,Cfc,, 



h GS 



-E N - 1 -e F 



h GS 



-E N + 1 +e F 



^E 

S' 



-e F 



h GS 



E N- 



-e F 



~ h GS 



(9) 



where, due to spin conservation, £ has the same direction than a and £' the same than a'. The first term is simply 
a constant. The second term (involving cl, _.Cfe )0 .) represents a potential scattering which does not involve any spin 
flip. These two terms are identical, to the ones appearing when building up an sd Hamiltonian from the Anderson 
Hamiltonian in the case of N = lE3. As it is usually done in that case, one can forget about these two terms which 
do not contain anything important for the physics we want to address. 

The interesting physics is included in the third term of (^|) which is a Kondo Hamiltonian 



H K - Efe,fc 

with exchange couplings 
J 



A+.-fi-A+.j; 



- e F 



^GS 



+ 



E N- 



- £F 



^GS 



A 2 



e F 



E N 
^GS 



E N - 1 +e F 



E N ' 

^GS 



(10) 



(11) 



3 



Hk is a spin-flip scattering Hamiltonian in which the 
GS of the scatterer flips its spin by means of X^^ and 
X^ (as shown in Fig. ||) and, at the same time, changes 
its total angular momentum. Both the difference between 
Mfl. and M^, and the correlation effects included in the 
tunneling spectral amplitudes, are the reason why one 
must use spin-flip Hubbard operators instead of the usual 
spin-ladder operators 



/ 



(M,S Z ) 



\ / 



|GS_ Z > 



(M',S z+ 2a) 



X 



|GS Z > 



FIG. 3. Simplest diagram representing the spin-flip scat- 
tering described by the Xe,-e = \GS^) {GS^ S \ Hubbard op- 
erators. The base line depicts the evolution of the momentum 
and spin of the scatterer (the QD). Dashed lines depict car- 
riers with spin a = ±1/2 at the leads. Spin conservation 
requires that a has the same direction than E. 

A crucial question is the sign of the exchange couplings 
(|TT|). With our definitions (Q), all the tunneling spectral 
amplitudes are positive. Therefore, the signs of the ex- 
change couplings are determined by the denominators. 
As, in the considered situation, the lowest energy corre- 
sponds to having N electrons in the QD, all the inter- 
mediate states 1 7), with N ± 1 electrons in the QD and 
=Fl electron at the Fermi level of the leads, have higher 
energy. Therefore, we have a very important result: TTr- 
has positive effec tive exchange couplings. So, a Kondo 
Hamiltonian ( |l0|) with antiferromagnetic couplings (|ll| ) 
has been obtained. The conclusion is: for any values 
of both N and spin, the QD in the regime 2 > v > 1 
presents Kondo physics. 



IV. EXCHANGE COUPLINGS AND KONDO 
TEMPERATURE 

Since the Hamiltonian ([l0|) is formally equivalent to 
the standard Kondo Hamiltonian for N = 1, the tem- 
perature dependence of the conductance (as any other 
measurable property) presents now characteristics simi- 
lar to the one found in those experimentsBu which have 
been interpreted in terms of only one electron. The im- 
portant issue to be discussed is the characteristic energy 
scale, Tk, which is determined by the antiferromagnetic 
couplings. 



J and Js depend on both the energy difference Eq S — 
E N±1 =p ef, and tunneling spectral amplitudes A_ ; £. 
The former is practically independent on both N and 
E. Therefore, it does not imply any significant difference 
with respect to the well known case of N = 1. However 
interesting physical differences appear due to A-^. Let 
us analyze the different situations in which the QD can 
be found: 

i) (Upper part of Fig. |l|) When the Zeeman interaction 
is large (for instance having a large inplane component of 
the magnetic field), the GS's are always compact states 



\GS») 



\C 



K ) 

N — KI 



\GS$) 



\n K + l 



N-K-H 



This is the simplest case in which 



X 



-1) 



d N _ K _ifd,K,i 



(12) 



(13) 



and similarly for the other Hubbard operators. We have 
also assumed, |7V— 1) = \C%_ K _ X ) and \N+l) = \C%+ K ) 
(in other case the quantum numbers of the gT operators 
must be changed accordingly) and the signs are taken ac- 
cording with definitions (g) . Due to the lack of correla- 
tion effects in the GS's, the tunneling spectral amplitudes 



A_ 



V N -k-i ; A_ 



•u- 



V K . (14) 



There is a difference between the antiferromagnetic cou- 
plings J, Jfl. and Jjj, due to the fact that N — K — 1 > K . 
It is easier to tunnel from the leads to spin up states 
(to = N — K — 1) within the QD because they are in 
the outer region of the QD while the first available spin 
down state (to = K) is in the inner region of the QD. 
This tunneling amplitude effect is-i^pf described by pre- 
vious models of Kondo at finite BEJ~t3, but has long been 
more broadly recognized both experimentally and theo- 
retically for QD with 2 > v > 1. 

Although in this regime the QD does not present any 
correlation effects, the antiferromagnetic couplings are 
still a function of total energies and tunneling amplitudes. 
Siiice the Kondo temperature depends exponentially on 
JB,it is a very sensitive magnitude with respect to many 
parameters as u>, u>q, B 7 V m etc. However, this is not dif- 
ferent from the caae_of just one electron in the QD in 
which experimentsEJ cl show Tk to be in the range of IK. 
In any case, one can predict variations in Tk when the 
regime changes as discussed below. 

ii) (Lower part of Fig. When the Zeeman interac- 
tion is small enough, therGS's are skyrmion-like states of 
topological charge one&3't3: 



\GS$) = \SK^±) 



(15) 



The Hubbard operators Xx;,sy have complicated, but an- 
alytical, expressions in terms of A^. 

When N is very large, correlation effects provoke that 
A_ £ tends to zero (for instance A_,s oc [NlnN]^ 1 / 2 for 
P = 0) implying a quenching of the antiferromagnetic 
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couplings due to the orthogonalization catastrophe. As a 
consequence, Kondo effect should not be observed when 
the number of electrons within the QD is large because 
the Kondo temperature is extremely small in this case. 

Even more interesting is the effect for a reduced num- 
ber of electrons, for instance N = 8 as both in Fig. [|, 
^ and in the experimentEll. In this case, correlation ef- 
fects do not destroy Kondo effect but reduce A_,£ up 
to a factor of twoED. The reduction of J and Js implies 
that Kondo temperatures are significantly smaller than 
for the case i) of compact states. In practice, one can 
move from the skyrmions regime ii) (lower part of Fig 
0) to the compact regime i ) (upper part of Fig [l]) by in- 
creasing an inplane component of the magnetic field. As 
a consequence of the above analysis, one should detect a 
clearly measurable increase of Tr- during this process. 

Hi) There is a rather unusual case, in which, for inter- 
mediate Zccman interaction and only for a few values of 
the magnetic field, the degenerate GS's are one skyrmion 
and one compact state, i. e.: 

\GS$) = \SK% iKt+ ) ; \GS$) = \C«I K+1 ). (16) 

This happens, for instance, in the last step to the right 
in Fig. H In this case, the two GS have spins differing in 
more than 1. Therefore, they can not be obtained from 
the same \N —1) by creating one electron with spin cither 
up or down. In other words, one of the two GS's has a 
tunneling spectral amplitude equal to zero so that, for 
tunneling effects there is only a non degenerate GS and 
Kondo effect does not occur. 

A practical difference between the compact and 
skyrmion cases resides in the probability of finding two 
degenerate (or almost) GS's. As shown in Fig. [I], 
the crossing regions are much closer to each other for 
skyrmions (lower part of Fig. [l]) than for compact states 
(upper Fig. |f|). To work with a magnetic field in which 
Kondo effect exists is more probable in the skyrmion 
regime than in the compact one. However, once Kondo 
regime is found in the compact regime, the effect is 
more robust because its Tk is higher than the one for 
skyrmions. 



V. CHESS BOARD BEHAVIOR OF THE 
CONDUCTANCE 



A very characteristic feature of some experiment: 




is an alternating high-low conductance sequence as a 
function of B for a given temperature and number of 
electrons. When N is varied in ±1, the high and low 
conductance valleys are interchanged. Since a common 
way of representing the experimental results is to use a 
color-intensity scale for the conductance as a function of 
both B and a gate potential .which varies N, the data 
present a chess board aspectUO. This occurs in a broad 
range of filling factors and number of electrons. 



In order to understand such an experimental behavior, 
let us start by applying our description for thej^jgime 
1 < v < 2. Since the chess board is observedE3E3 in 
samples with large w, skyrmion- like states are high en- 
ergy excitations. GS's are restricted to be given just by 
one configuration, i. e. [S^'Y\C^_ JC ) states, where j 
is an integer (for practical purposes one can concentrate 
in the cases j — and I). At first sight, the considera- 
tion of both \C§_ K ) and S^\C§_ 

k) seems unnecessary 
because these states have an energy difference equal to 
the Zeeman coupling so that they are never degenerate. 
However, in a broad range of magnetic fields, these two 
states are precisely the two lowest energy states (see e.g. 
Fig. §. 



q(-)ip k+1 ^ 

5> N-k-1> 




Conductance 




B 

FIG. 4. Schematic evolution with the magnetic field of the 
two lowest energy states of a QD (upper part) and the con- 
ductance (lower part) in the range Tk > T > AE 2 /Tk (see 
text). 

In the two level system approachEll in which the two 
lowest states are separated by an energy splitting AE, 
Kondo-like behavior appears only when the experimen- 
tal temperature is between (AE) 2 /Tk and Tk (provided 
AE < Tk) in order to have significant occupation of the 
two states. The quantitative application of the model 
would require the computation of both the energy split- 
ting AE and Tk which, as mentioned above, are very 
sensitive to many experimental parameters. Instead, a 
general understanding of the chess board is obtained from 
the following qualitative explanation depicted schemati- 
cally in Fig. 0: when magnetic field is varied, the two 
lowest levels are in two alternating situations depend- 
ing whether the system is around or far from a cross- 
ing of GS's. Around a crossing, the two lowest levels 
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correspond to states ICj^^-i) and \C^_ K ) while far 
from the crossing they are \C$_ K ) and S^\C^_ K ). In 
the former case, the two lowest levels can be connected 
by the tunneling Hamiltonian through low energy states 
|JV±1) with A±,s = VkVn-k-i- In this case, the system 
presents high Kondo-like conductance. In the second case 
of being far from the crossings, the two lowest states are 
\C§_ K ) and S { - ) \C§_ K ). The coupling J between the 
two states is very weak for two reasons: first, the connect- 
ing states |JV± 1) are highly excited states implying very 
large denominators in J and, second, since S^' \C^_ K ) 
is a linear combination of different configurations, A± ( e 
is significantly reduced by factors proportional to 1/N. 
These two effects produce, in this region, a very small J 
and, consequently, an exponentially negligible Tk so that 
the measured conductance is not Kondo-like but instead 
is quenched. This originates the alternating behavior ex- 
perimentally observed for fix gate voltage (i. e. num- 
ber of electrons in the QD) and varying magnetic field. 
The chess board aspect also implies alternating low-high 
conductance regions for fix B and varying gate voltage, 
which is due to the fact that crossings for N ± 1 electrons 
occur for magnetic fields roughly midway from crossings 
for N electronaiZI. 

Since Kondo-like behavior only occurs when the exper- 
imental temperature is in the range between (AE) 2 /Tk 
and Tk , a clear prediction of our scheme is that the highly 
conducting regions of the chess board would become nar- 
rower for decreasing temperature due to the lower limit 
condition. 

The sequence of alternating high-low conductance has 
been observed also at very low magnetic-field (B ~ 0.5T) 
and large number (TV ~ 50) of electronsEl So, a complete 
understanding of the experimental situation requires the 
extension of our framework to v > 2. We are currently 
involved in this task which implies to take into account 
GS's more complicated than simple compact states in the 
lowest Landau level. 



VI. SUMMARY 

We have presented a general description of the Kondo 
effect for any number of electrons in a QD at filling factor 
1 < v < 2. Collective spin effects of the N electrons are 
described in terms of a set of spin-flip Hubbard operators 
which also change the third component of the total an- 
gular momentum. Our conclusion is that, for any num- 
ber of electrons within the QD, the spin- flip scattering 
of carriers only produces transitions between two (either 
compact or skyrmion-like) states of the QD with spins 
differing in 1. This process can be described by a Kondo 
Hamiltonian with exchange couplings that, we show, are 
antiferromagnetic and depend on both correlation and 
tunneling amplitude effects. As a consequence, the sys- 
tem presents a Kondo behavior with an experimentally 
accessible Tk ■ The increase of Tk with an inplane com- 



ponent of the magnetic field should allow the detection of 
the transition from a skyrmion-like regime to a compact 
state regime. Finally, we present a qualitative explana- 
tion for the chess board aspect of the experimental con- 
ductance when represented in a grey scale as a function 
of both the magnetic field and the gate potential affecting 
the quantum dot. 
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